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Abstract 

The physical aspect of a general perturbation theory is explored. 
Its role as a physical principle for understanding the interaction among 
matter with different levels of hierarchy is appreciated. It is shown 
that the generic perturbation theory can not only be used for under- 
standing various electronic phenomena including the nature of chem- 
ical bonds but also serve as a unified theme for developing general 
electronic structure theories and calculation schemes. In particular, 
a standard electron correlation approach is suggested and established 
according to this law. 
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Perturbation theory is regarded as one of two major approaches for ap- 
proximately solving quantum many-body problems. However, its deeper 
physical aspect is far more than it is at present being used just as a math- 
ematical apparatus for solving the complicated issues. All the fundamental 
laws in physics are variational in nature, including the Schrodinger equa- 
tion in quantum mechanics. Nevertheless, the perturbation theory provides 
a basic principle that governs how matter with different levels of hierarchy 
interacts. In fact, a general perturbation method itself contains two ingre- 
dients. On one hand, the degenerate or near-degenerate situation is not a 
perturbation at all but actually constitutes a strong physical interaction. On 
the other hand, the non-degenerate case is a real perturbation in the common 
sense. We believe that, it is this physical mixing with equal or near energies 
that governs the interaction among matter with different levels of hierarchy. 
Of course, it is also an elementary physical law based on which a unified 
chemical bond theory can be built. 

Electrons are quantum mechanical entities which possess wave-particle 
duality. The binding process of the electrons associated with some atoms, 
or equivalently, the interaction of atomic orbitals for the formation of a 
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molecule, can be regarded as a wave interference phenomenon. The interac- 
tion of intra-atomic orbitals with the same or near energies is the Pauling's 
hybridization process, which determines the spatial orientation of chemical 
bonds, while the interaction of inter-atomic orbitals with the same energy or 
near energies decides the actual formation of chemical bonds. These are the 
nature of chemical bonds PUI2]- The immediate benefit for recognizing this 
near energy principle in the determination of the chemical bonds is that it 
gives a better understanding of many very important structural concepts such 
as multi-center chemical bonds, multiple chemical bonds, resonance struc- 
ture, Walsh diagrams, and avoided crossing. It can incorporate these differ- 
ent concepts into single qualitative theoretical framework [H El 13 El E] ■ 
Furthermore, in addition to serving as a fundamental physical law for 
understanding how matter interacts, the generic perturbation approach can 
provide a unique and most powerful mathematical device for quantitatively 
investigating the electronic structure of molecules including molecular mate- 
rials and biomolecules 0. We are going to have a harmony between theory 
and computation. 

Energy scale principle in the Rayleigh — Ritz variational approach 
The Rayleigh- Ritz variational method is most commonly used for solving 
eigenvalue problem in quantum mechanics. Its relation to the general pertur- 
bation theory has been the subject of analysis for many decades, which in- 
cludes the study of the complicated issues concerning the perturbation theory 
for the linear operators [3 EH E] ■ However, the utilization of this relation, 
or its role as a guidance, in establishing the electronic structure calculation 
schemes has not been fully explored or completed yet. A related universal 
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formalism for the quantum many-body systems is still lacking, which is of 
paramount importance for our investigation ■ At a first glance, as long as 
the reference Hamiltonian which produces the basis functions is made as close 
as possible to the full Hamiltonian, then the dimension for the Rayleigh-Ritz 
variational expansion will be made as small as possible. In addition, if the 
basis functions have the closest energies of the reference Hamiltonian, then 
they will have the strongest mixing and make the greatest contribution to 
the combined states; while the others with larger energy differences will have 
smaller or even negligible contributions. These are the situations we qualita- 
tively discussed above for the elementary perturbation theory. We term this 
as the energy scale principle in the Rayleigh-Ritz variational approach, 
(a) Molecular fragmentation and basis set construction 
The basis set approach is a most popular and natural way for solving the 
single particle equations such as the Hartree-Fock equation for a molecular 
system. Physically, it reflects a composite relation between the molecule and 
its constituent atoms. To have an overall accurate electronic structure calcu- 
lation, the first necessary step is to get the reliable and converged molecular 
orbitals [T3] . 

However, since the current basis functions like contracted Gaussians which 
are most commonly used are primarily a reflection of electrons in single atoms 
in the molecule, it leaves the perturbed part of the molecular Fock operator 
very large. That is why the polarization functions, including some expanded 
ones such as the correlation consistent basis sets, have to be introduced to get 
good computation results [T4*l ll5|. Nevertheless, the 0(N 4 ) scaling, where N 
is the number of basis functions, has become a major bottleneck in quantum 
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chemistry calculation, especially for the large systems. 

To overcome this difficulty, the energy-scale principle described above 
can be helpful. If we construct the basis functions which are the reflection 
of molecular fragments so that the corresponding reference Hamiltonian is 
as close as possible to the whole molecular Fock operator, then the dimen- 
sion of basis set expansion can be made as small as possible. This is going 
to be a challenge work but will be mathematical in nature. The basis set 
superposition effects (BSSE) is an example [To] . 

A similar situation occurs in the quantum molecular scattering calcu- 
lation, where the channels are used as the basis functions for solving the 
Schrodinger equation, or its integral form, Lippmann-Schwinger equation 
with proper boundary conditions. Since there are often very large differences 
between the channels and the scattering waves for the whole reactive system 
in the interaction regions, the dimension for their expansion is particularly 
large. This causes the quantum scattering calculation to be prohibitively 
expensive for all but the smallest systems. The ideas from the perturbation 
theory can obviously be utilized for remedying this deficiency [T7] . 

(b) General multireference electronic structure theory 

To get a final accurate solution to the Schrodinger equation for the many- 
electron system with a non-separable two-body Coulomb potential, it is most 
likely that we have to go beyond the single particle description [THl EH I2D1 
I2T] . Mathematically, the full configuration interaction (FCI) gives exact 
answers [22] ■ However, it is computationally prohibitive and possibly will 
never be strictly realized. The energy scale principle described above can 
also be applied in this configuration level. A general electronic structure 
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theory must be multiconfiguratiorial or multireferential in nature |2H1 
1213 123 I2H I2H E01- First, there exists a strong configuration mixing, for 
example, at transition states, for excited states, and for multiple chemical 
bonds. The concept of exciton introduced in solid state physics also belongs 
to this case [UJET1E2]- Second, the degenerate configurations are often the 
case for the stable open-shell systems. Third, if we want to treat the ground 
state and the excited states simultaneously, we have to include the corre- 
sponding reference states in the same model space. Finally, the separation 
of correlation into the static and dynamic parts, which corresponds to the 
near degenerate and the perturbed situations, really has chemical structure 
signature. Therefore, among all the correlation approaches developed so far 
for the electronic structure, the type of multiconfigurational self-consistent 
fields (MCSCF) with perturbation or coupled-cluster expansion corrections 
should be the most appropriate one and works in the right direction. To solve 
the issues such as proper selection of the configurations for the model space, 
the efficient treatment of the dynamic correlation, and the avoidance of the 
intruder states, we need not only a mastery of current quantum many-body 
theories but also their further development pffi E3] . 

The density functional theory (DFT) is one of the most powerful compu- 
tational schemes in the study of the electronic structures for the molecules 
and solids jHSHSHl- However, the importance and necessity of the separation 
of the correlation into a static part and a dynamic one is also indicated in its 
calculation of the highly charged ions and in its treatment of the transition 
states for reactions [HZ1EE]- Even though the time- dependent DFT has been 
developed to address the excited state issues, it seems that the overall den- 
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sity functional theories are still within the framework of quantum mechanics 
based on the state vectors and operators [SHI ■ The density is one of the most 
fundamental physical quantities which is used for specifying a system. 

One most generic quantum many-body approach has been suggested and 
established within the perturbation theory ^3] • Since the molecular orbitals 
define a reference Hamiltonian for the interacting many-electron system jlUl , 
it is crucial to choose the appropriate ones to make the reference Hamilto- 
nian as close as possible to the whole Hamiltonian. In this way, the compu- 
tational burden relevant to the correlations will be alleviated. Nevertheless, 
the correlated molecular orbitals, determined from the Dyson equation, pro- 
vide this candidacy. It is a most general energy eignequation for generating 
the molecular orbitals, covering the Hartree-Fock and Kohn-Sham equations 
as special cases. In practice, it gives the best smg/e-particle properties such 
as electron affinities and electron ionization potentials. Meanwhile, the one- 
particle Green function method can go beyond the single-determinant level 
and study the iV-body properties such as the transition states and total ener- 
gies including excitation energies related to the configuration mixing. Since 
the self-energy operator can be approached in a systematic way, the multiref- 
erence perturbation theory, combined with the single-particle Green function 
formalism, will henceforth furnish a most powerful approach for studying the 
static and nonstatic correlations of the interacting electrons [T3|. 

(c) General pseudopotential theory 

The concepts of effective core potentials (ECP), pseudopotentials, or 
model potentials (MP) are those of the most significant developments in the 
fields of the electronic structure for molecular and solid state systems [JTJ 
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W2\ . They treat valence electrons only, leaving the core electrons and the 
nucleus as a whole charge entity and therefore reducing the number of elec- 
trons as well as the corresponding overall size of the basis sets employed for 
the computation. It is important when we study the electronic structure for 
the large molecules or the inorganic molecules containing heavy elements. A 
type of pseudopotentials most commonly used for the solid state calculations 
is the so-called norm conserving pseudopotentials [42J. In addition to having 
the same energies for the valence states, their pseudo valence wavefunctions 
are equivalent to the wavefunctions for the valence states obtained from the 
full electron calculations outside a cutoff radius. The pseudopotentials con- 
structed in this manner share the same scattering properties as those of the 
full potentials over the energy range of the valence states. The practical 
implementation of the various pseudopotentials has also demonstrated the 
importance of choosing a correct size of the core or the range of the va- 
lence electrons for the accurate pseudopotential computation in order that 
the core-valence correlation or core polarization can be neglected. Obviously, 
the physics behind this separation of the valence and core states is the energy 
scale principle we described above applied in the level of atomic orbitals. 
After realizing this principle, however, we might establish a more general 
pseudopotential theory. We are planning to do so within the framework of a 
perturbation theory so that the most flexible and accurate pseudopotentials 
or effective core potentials can be obtained. They can be used in different 
chemical environments and work for both ground and excited state problems. 
The final goal is to make the effective core potentials to be a routine for the 
calculation of the electronic structures for the large molecules, the inorganic 
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molecules containing heavy elements, and the solid state systems. 

(d) Molecular fragmentation and combined quantum mechanical and 
molecular mechanical ( QM / MM) approach for electronic structure of 
large molecules 

A combined QM/MM approach has become very popular in recent years 
in the investigation of, for example, the chemical reactions in solutions and 
in enzymes |44j . The basic consideration is that treating a full collection of 
electrons for the whole system explicitly is not only unrealistic but also un- 
necessary. In the first place, the electronic charge redistribution induced by 
a chemical reaction is very often limited to a small region due to the length 
scale issues such as a finite range of interaction or natural charge distribu- 
tion. Second, the quantum exchange effect for the electrons is finite range, 
and there is no exchange interaction among the electrons with a long dis- 
tance |43[ IIHI El EH ES]. This permits a partition of the whole system into 
an active part and an inactive one without any charge redistribution. The 
former has to be described quantum mechanically since it possibly involves 
bond breaking and making, while the latter can be described by molecular 
mechanics because it merely serves as a classical electrostatic environment 
for the active site [SHI- This combined QM/MM description has shown re- 
markable successes in studying the electronic structure and reactivity of large 
molecules. However, challenges remain. One of the major obstacles for the 
applications is in the proper treatment of the boundary region where the 
cut has to be for a covalent bond. Currently, there are two approaches to 
this problem. The one introducing link atoms along the boundary is severely 
limited and cannot be applied to treat a large variety of different chemical sys- 
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terns. In addition, it artificially brings additional forces into the system and 
therefore complicates the problem. The other kind like local self-consistent 
field methods seems reasonable but it is still more empirical. In order to uti- 
lize this kind of combined QM/MM methods for investigating the electronic 
structure and molecular dynamics in a larger domain of fields, we need to 
develop a more generic ab initio approach. We believe that the energy scale 
principle discussed above can play a key role here. It is not only the prin- 
ciple according to which the atomic orbitals including valence ones interact 
along the boundary but also the law based on which a systematic approach 
for constructing the correct charge distribution or the force fields along the 
boundary can be established. This is also the key for a more sophisticated 
or finer treatment of the quantum region including its electron correlation. 

In summary, the energy scale principle for the hierarchy of interacting 
matter is identified. Not only can it be utilized as a universal law for under- 
standing how matter interacts at different levels but also the relevant pertur- 
bation formalisms, including the ones pertaining to the Rayleigh-Ritz varia- 
tional expansion, can serve as the foundation for building various and most 
powerful electronic structure calculation schemes even for the large molecu- 
lar systems [U]. In particular, a standard electron correlation approach or a 
quantum many-particle theory in general have thereby been established [T3] . 
Obviously, they can also be employed to develop a generic theory for the 
intermolecular forces so that the important issues such as the interplay be- 
tween chemical bondings and intermolecular forces can be investigated [5T] 

E2|. 
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